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^^ ■ PACS. 05. 45. Ac - Low-dimensional chaos. 

.^.^ ' PACS. 51.20.+d - Viscosity, diffusion, and thermal conductivity. 

\^ ' Abstract. - A consistent description of a shear flow, the accompanied viscous heating, and 

the associated entropy balance is given in the framework of a deterministic dynamical system, 
where a multibaker dynamics drives two fields: the velocity and the temperature distributions. 
In an appropriate macroscopic limit their transport equations go over into the Navier-Stokes 
and the heat conduction equation of viscous flows. The inclusion of an artificial heat sink can 
stabilize steady states with constant temperatures. It mimics a thermostating algorithm used 
, in non-equilibrium molecular-dynamics simulations. 

cn 

o 

f^ , Introduction. - In the last years, there has been an increasing interest in modeling 

^^ ' transport phenomena by low-dimensional, deterministic dynamical systems p|-p^ . Multibaker 

^p I maps p|-[l3| appeared to be the simplest models of this approach. They provide an opportunity 

to derive the equations of non-equilibrium thermodynamics from an underlying dynamics 

without using the concept of particles. The strongly-chaotic mixing properties of these two- 

G ^ dimensional maps seem to be sufficient to ensure consistency with the entropy balance equation 

^ I of thermodynamics provided that a properly chosen coarse-grained entropy and a macroscopic 

r '~j ' limit are taken |Q,|l^,Q . Previous work successfully described the phenomena of diffusion |Q,g| , 

rN I conduction in an external field f^, 0,|^, chemical reactions g], thermal conduction ||lO|, and 

j^ ' cross effects due to the simultaneous presence of an external field and heat conduction |ll|] by 

means of multibaker maps. Not only stationary, but also transient states could be addressed 

[^, 11, 1^. Here we add to this list the phenomenon of shear flows and the accompanying 



viscous heating. The interest of this is to clarify how the shear rate enters the expression 
for the irreversible entropy production. After all, by the definition of local equilibrium the 
macroscopic fiow profile does not appear in the entropy balance. Thermodynamic averages 
contain only deviations from the average streaming velocity. 
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Fig. 1 - The multibaker map. a) The chain of iV + 2 baker cells of height a and of unit width 
modeling an intersection of the fluid along the x axis. The arrows indicate the average velocity Vm of 
the fluid at positions x G [am, a{m + 1)]. b) The 'microscopic' dynamics of cell m (0 < jn < A*' + 1). 
All notations are explained in the text. 



We consider a viscous fluid driven by applying a shear forcing to two parallel walls of 
separation L {cf. (l4[^). The fluid is assumed to be incompressible. The upper wall moves 
in the y-direction with a constant velocity U (flg. ||a). We restrict ourselves to a discussion of 
(possibly non-stationary) laminar flow. Due to translation invariance the velocity v is always 
directed in the y direction, and it only depends on the x coordinate. Similarly, the temperature 
T is independent of the y and z coordinates. There is no external pressure gradient applied, 
and the internal pressure forces due to eventual temperature changes are expected to be 
negligible. For sufficiently long times the well-known linear velocity profile v* (x) = Ux/L is 
approached as asymptotic steady state. The spatial behavior of the temperature distribution 
T{x, t) depends on the boundary conditions. Initially T is typically increasing in time due 
to the viscous heating. For an adiabatically closed fluid, the asymptotic distribution T* is 
spatially constant and grows linearly with time. Inclusion of a heat sink can remove this time 
dependence such that the asymptotic state becomes stationary. In a hydrodynamic setting 
the sink acts only at the boundaries of the system, but an artiflcial, spatially-uniform sink 
can also be implemented. It mimics the so-called SLLOD ^^ algorithm of non-equilibrium 
molecular dynamics. 

A multibaker m,ap for shear flow. - We model an intersection of the fluid along the x 
axis by a multibaker map. It consists of A^ -f 2 cells of size unity times a (fig. |l|a). The 
cells m = and N + 1 represent the boundary, and cells m = 1, ..., iV, the bulk of the fiuid, 
whose width is Na = L. In each cell there is a momentum- like variable p defined besides 
the position. The multibaker dynamics advects two fields; the 'microscopic' velocity v{x,p;t) 
and the 'microscopic' temperature T{x,p;t). The appropriate (c/. below) averages Vm and T^ 
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over cell m are called the coarse-grained fields. In the macroscopic limit they go over into the 
hydrodynamic velocity field and the thermodynamic temperature, respectively. 

The two-dimensional multibaker dynamics is defined as follows (fig. ^b). After each time 
unit T every cell m = 1, • ■ • A^ is divided into three bands of heights ag, as and ag, such that 
g + s -\- g = 1. The outermost ones are mapped onto a column of height a and width g in 
cells m + 1 and to — 1, respectively. The middle one stays in cell m, where it is transformed 
into a column of width s. In all cases the area is preserved such that an initially constant 
phase-space density g remains constant in time, reflecting the incompressibility of the fluid. 
Hence, the density can be identifled with the density of the fluid. 

The velocity and the temperature fields are advected by the multibaker dynamics. In 
contrast to the density g they typically evolve in time. Starting with a coarse-grained velocity 
distribution Vm along the chain, the updated values w^ after time r become 

v'jn = Vm + g{Vrn-l + ^m+l - 2i;„i), (1) 

This update expresses momentum conservation in the x-direction. A portion s = I — 2g of the 
original momentum (velocity) remains in cell m, and portions g of the momenta (velocities) 
of the neighboring cells fiow in. 

The temperature equation follows from the energy balance. The full energy Cm of cell m 
after coarse graining is proportional to the sum of the temperature T„i and the translational 
kinetic energy v'^/2, i.e., 

where C is a constant. The update of energy is due to an in- and outflow of energy from the 
neighbors 

^''"^' ' CT:,, = {l-2g)(^-^ + CTm]+g(^-^ + CT,„^+g(^-k^+CT,^+^. (3) 



2 " ' "' \ 2 

The action of a thermostat can, however, lead to a change of this kinetic energy. This is 
modeled by introducing a local source term q^ , and multiplying the energy terms by a factor 
[1 -I- Tqm], i-c, by setting e^ — > e^ [1 -|- rgm] after every update (^. Rearranging eq. (H) 
subjected to this additional factor, and using relation (ll]) one obtains 

TL = \ T„i + g(rm-i + T,n+1 - 2T„i) + ^[(Wm-1 ~ v^f + («„+! - V„if] \ [I + Tq^], (4) 

which can also be written in the form of the balance equation 

rpl rp •(«') _ Aw) 

^m ™ — n _ -^m+l Jm , ■. 



T 

with the discrete "heat" current 



„2 7^ _ 7^ 

Ji^ = — 9 : ■ (6a) 



and the source term 



1m rjil 



l + rqm ™ 2tC 



(Wm-l - Urn) + [Vm+l - VmY . (6b) 
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The first contribution to the source reflects the action of thermostating, and the second one 
the effect of viscous heating of the fluid. 

In the absence of coarse graining, a repeated application of the chaotic, mixing niultibaker 
dynamics leads to a fractal distribution of the temperature and the velocity fields. Similar to 
the treatment of diffusive mass and heat transport M, ni| the emergence and proliferation of 
these structures lies at the heart of a consistent dynamical-systems treatment of static and 
transient transport phenomena. 

Entropy and entropy balance. - The Gibbs entropy of a niultibaker system is defined as 
the information-theoretic entropy, i.e., as the phase-space average of lii-^g/ g*). The reference 
density g* in the single-particle phase space of the multibaker map is expected to depend on 
the temperature. We take the choice g* — gT^ , where 7 is a constant ex pon ent. Since the 
density g is constant in space and time the Gibbs energy of cell m becomesPj 



S'(G) ^ 



gi 



/dxdp \nT{x,p) (7) 



The coarse-grained entropy is defined in an analogous way as 

Sm^ajglnTm- (8) 

It is based on the cell-averaged, i.e., coarse-grained value of the temperature T{x,p). As 
announced, the average (streaming) velocities Vm do not enter the definition of the entropy. 

For the purpose of deriving the time evolution of entropies, it is useful to choose and initial 
condition with uniform densities in every cell. The coarse-grained and the Gibbs entropy then 
initially coincide. After one time step, however, the Gibbs entropy changes due to the fact 
that the T field takes different values in the neighboring cells, viz. 

5,^f )' = jga {g ln[T,„_i (1 + rq,-,,)] + (1 - 2g) ln[T™ (1 + Tq„,)] + g ln[T„,_i (1 + rg™)]} 

= gia \ ln[r„ (1 + Tg„,)] + 5 In ^^ - 3 In — ^^ \ . (9) 

t Jm -Lm+1 J 

On the other hand, the coarse-grained entropy after one time step is 

S'„, = ~a'jglnT:^. (10) 

Since it only depends on averages in small volumes in the configuration space, the coarse- 
grained entropy is considered as the analogue of the thermodynamic entropy M, ni|. Its 
temporal change can be decomposed as 



A Q Q' — Q c(G) _ c{G) / q i _ q(G) \ _ ( n _ c{G)\ 

^-^^m ^m ^^n Jm Jm \^m Jm } \Jm Jm } 

ar ar ar ar 

and information-theoretic arguments |7,11| lead one to identify 

A <? qid) _ q(G) 



ar ar 



(11) 



(12a) 



iXibm ('-'m i'm ) ('-'m »^r?j ) H ^V.^ 



ar 



( )The Boltzmann constant ks and a constant additive constant are suppressed here; c/. |13] for details. 
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with the entropy flux and the rate of entropy production, respectively. Note that the second 



term of the numerator of eq. (121:) vanishes due to the choice of initial conditions. 

Inserting eqs. (H) and (p^) into eq. (121;), the rate of entropy production AiSm/{aT) per 



unit volume and time is found to be 

JQT 



A"? 



ln{ ^:^{1 + Tq^) -.gin— .gin— ^- 



(13) 



It does not depend on the source term g,„, but only on the values of the field T in cell m 
and its neighbors. Entropy production arises in this model from (i) an explicit evolution 
of the macroscopic temperature profile, and (ii) from mixing of regions with different local 
temperatures. The emergence of self-similar structure in the temperature distribution T{x,p) 
can lead to a non-vanishing steady-state entropy production. Note that on the level of this 
discrete relations the entropy production does not yet contain the velocity distribution Vm- 
The effect of shear flow is only implicitly arising from the splitting of the full kinetic energy into 
a translational and an irregular part that depends on the implicit choice of a spatial resolution 
when writing eq. (Q). In the Gibbs entropy the local velocity v{x,p) enters implicitly through 
the definition (g) of the temperature, while for the coarse-grained entropy only the averages 
on the cells {i.e., w,„) can enter. 

The entropy flux lS.f^Sm/ {o-t) can be written as the sum of the discrete divergence of the 
entropy current 

£«^ = -^7Pln%±i (14a) 



and the flux 



into the thermostat. 



^'i^^ = IQqm (14b) 



The macroscopic limit. - Establishing a discrete momentum, energy and entropy balance 
is not sufficient to motivate the thermodynamic relevance of a dynamical systems model 
transport. As argued in H a full consistency can only be found in a continuum scaling 
limit (the m,acroscopic lim,it) where the time evolution equations have to coincide with the 
relations known from irreversible thermodynamics, irrespective of the detailed prescription 
of how to choose the discrete time and space units needed to define the local equilibrium. 
The macroscopic limit corresponds to a <^ L, TV ^ 1, and r is much smaller than typical 
macroscopic time scales (for instance the viscous relaxation time). Formally it is defined as 
a, T — > such that the spatial coordinate x — am is finite. Taking this limit will be indicated 
by the arrow — >. 

In order to find the macroscopic limit dtv — v d^v. of the coarse-grained velocity dynamics 
eq. (El) it is required that go?' /t takes the finite value v in the limit, i.e., 



a" 



(15) 



For a vanishing pressure gradient eq. (|l|) leads to the Navier-Stokes equation for the laminar 
shear flow [[l4[ , where v is the kinem,atic viscosity. 

Considering the macroscopic limit of eq. (B) one obtains the temperature equation 



,2^ ^ 



dtT^:^d'^T+-{d,Vyr + Tq. (16) 
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For q ~ this is exactly of the type known from hydrodynamics |l4], where, however, the heat 
diffusion coefficient k appears in front of the second derivative. Hence, in the muhibaker map 
v also governs heat diffusion, i.e., n — v. In the hydrodynamic expression, the coefficient of 
the term expressing viscous heating is v /cv, where Cv is the specific heat at constant volume. 
Thus, we have to identify the constant C with the specific heat, as also expected from a 
physical interpretation of eq. (0) . 

Evaluating the macroscopic limit of eq. (Il^) one gets for the entropy density 5,„/a — > s 



^^^ dts = a(^") - Vj(^) + $(*'') (17) 



ar 

which coincides with the thermodynamic entropy balance [ [l4[ if $(*'') = in the bulk. For 
the irreversible entropy production we find 

^iSm (irr) ( OxT\ VQ 2 /•i o\ 

'^ ^=ive\-^\ -^^{p^v) . (18) 



a^ 



It is consistent with thermodynamics if v^q corresponds to the heat conductivity A of the 
flow. Since in general A = k.qCv^ and since k. = v yo. our case, we conclude that also 7 amounts 
to the specific heat in our model, i.e., c^ = C — 7. 

For the entropy current jm one finds in the macroscopic limit 

£^^-/^^(-) = -7-P^ = -A^. (19) 

Also this relation fully agrees with its thermodynamic counterpart W% . The heat conductivity 
A = KQCy = VQ^ appears in front of the logarithmic derivative dxT /T without possibility to 
adjust free parameters, thus demonstrating the full consistency of the results with irreversible 
thermodynamics. Note that eqs. (fTh^diy) are valid at any instant of time. 

Conclusion. - We have established a simple multibaker model that faithfully models 
hydrodynamic and thermal properties of shear flows. The model is based on an area-preserving 
multibaker dynamics. No phase-space contraction was needed. Thermostating manifests itself 
in the appearance of the a source term q that only influences the heat equation and the entropy 
flux. By a proper choice of q{x, t) one can stabilize any temperature profile T*{x) as a steady 
state. For instance the choice q = v{dxv) /{CT) eliminates viscous heating, and ensures a 
T* =constant steady state for a uniformly thermostated system. Alternatively, a parabolic 
stationary temperature profile is obtained when the source terms only act at the boundaries 
in order to prescribe a constant temperature at the two ends. In either case, the source term 
$("') leads to an entropy flux into a thermostat, and hence requires a generalization of the 
local relations of classical irreversible thermodynamics. 

The main interest of the model lies in the light it sheds on the origin of viscous heating in 
deterministic models of transport. As in the SLLOD algorithm it arises from the emergence 
of fractal structures. In contrast to that model where they are due to the joint action of a 
driving force and a Gaussian thermostat introduced by heuristic arguments into the equations 
of motion, the present model also generates the structures for more typical physical settings 
of transport driven from the boundaries. It identifies the structures as arising from the mixing 
of regions with different local temperatures that is exponentially proliferating to smaller and 
smaller scales for a driven system. In the present model one can explicitly follow the analog of 
an energy cascade in turbulence where the kinetic energy of a macroscopic flow is distributed 
to finer and finer scales until reaching the Kolmogorov scale where it has to be considered 
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as contributing to the non-directional motion and leads to viscous heating. It is exactly this 
mechanism that also leads to the appearance of the macroscopic shear rate in the expression 
of the irreversible entropy production. 
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